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Abstract 

We attempt to understand the final-state interaction in the two-body non- 
leptonic decay of a heavy particle for which many multibody {N > 3) decay 
channels are also open. No matter how many multibody channels couple to 
the two-body channels, the analyticity of the S-matrix relates the phase and 
the magnitude of the two-body decay amplitude through a dispersion relation. 
In general, however, the phase cannot be determined by strong interactions 
alone. The dispersion relation requires on a general ground that the final- 
state interaction phases be small for the two-body decay amplitudes when 
the initial particle is very heavy. We then analyze the final-state interac- 
tion phases in terms of the s-channel eigenstates of the S-matrix and obtain 
semiquantitative results applicable to the B decay with a random S-matrix 
hypothesis. We use the high-energy scattering data and the dual resonance 
model as a guide to the relevant aspects of strong interaction dynamics at 
long and intermediate distances. 

PACS numbers: 11.55.Fv, 11.80. Gw, 12.40.Nn, 13.20.Fc, 13.20. He 
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I. INTRODUCTION 



The final-state interaction in the nonleptonic weak decay is difficuh to estimate when 
a large number of multibody channels are open. While the short-distance final-state inter- 
action is small and its computation is noncontroversial, we have little understanding, 
theoretically or experimentally, of the long-distance final-state interactions. 

The long-distance final-state interaction phases were computed from the high-energy 
Regge exchange amplitudes in the elastic rescattering approximation . The experimental 
data appeared in favor of large final-state interaction phases at least for the D decay P|. 
It was asserted that the measured phases of the two-body D decay amplitudes can be 
reproduced in the elastic approximation to the final hadron interactions [^. However, it is 
fairly obvious from an analysis of the partial-wave unitarity with the diffractive scattering 
Q that the elastic approximation cannot be justified at high mass scales, for instance, in 
the B decay where the two-body final states can couple to a very large number of multibody 
final states. Even when we are interested in the final-state interaction phases of the two- 
body channels alone, we cannot determine them without knowing the coupling of the two- 
body channels to the multibody channels. The Regge amplitudes alone do not provide 
all necessary pieces of information. Actually, strong interactions and CP-conserving weak 
interactions are entangled in the decay phases. In this paper we shall make a modest attempt 
toward understanding of the inelastic final-state interactions. Because of the limitation in 
the numerical computation of the long-distance effects, we are able to present our results 
only in a semiquantitative way. 

We present two approaches here. The first one uses the analyticity of the S-matrix. For 
the two-body decay, the phase and the magnitude of amplitude are tightly related to each 
other by a dispersion relation no matter how many multibody channels couple to two-body 
channels. The same phase-amplitude relation does not hold for the multibody decay. So far 
no theorist has ever attempted to study the correlation between the phase and the magnitude 
from this aspect. No dynamical assumption is introduced in this approach. 

The origin of difficulty in the final-state interaction at high mass scales is in that so many 
channels are open and communicate with each other. In the second approach we analyze 
the decay in terms of the s-channel eigenstates of the S-matrix and treat a large number of 
open eigenchannels statistically by introducing a randomness hypothesis |]^ . While we give 
up much of numerical predictability in this approach, we are still able to see general trends 
in the final-state interaction at high mass scales. Both approaches lead us to conclude that 
the long-distance final-state interaction phases should be small for two-body heavy hadron 
decay. Though our conclusion favoring small final-state interactions may be in line of some 
of the existing literature, our method and picture are completely orthogonal to them. 

In Section II, after a brief review of the analyticity of the decay amplitude into general A^- 
body channels, we derive for the two-body decay a dispersion relation which relates the phase 
and the amplitude through the Omnes-Mushkelishvili integral Using this dispersion 
relation we separate from the physical decay amplitude the final-state interactions below an 
arbitrarily chosen timelike energy scale. We see in this form that a final-state interaction 
phase of any origin cannot persist to very high energies. 

In Section III, we shall study the final-state interaction phases from the viewpoint of 
the eigenphase shifts of the strong-interaction S-matrix. We write the hadron scattering 
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amplitude in terms of the eigenphase shifts, and classify the eigenchannels into the resonant 
and nonresonant ones according to the dual resonance model. Then we express the final- 
state interaction in terms of the eigenphase shifts. We study the high-energy behaviors of 
the eigenphase shifts from the scattering data. 

In Section IV we introduce the dynamical postulate that the composition of eigenchannels 
is statistically random when very many of them exist in degeneracy. We make a quantitative 
estimate of the decay amplitude phases within the limitation of the method. 

Finally in Section V, we apply our findings to the actual D and B decays. For the D 
decay, the elastic approximation combined with the Regge asymptotic behavior is not allowed 
in determining the final-state interaction phases. The random-phase method is probably a 
poor approximation. The observed large phase difference between different isospin channels 
can be accommodated but not predicted. The random phase approximation has the best 
chance in the B decay where the final hadron multiplicity is high. All two-body decay 
amplitudes are dominantly real up to CP violations in the B decay. If the color suppression 
exists prior to final-state interaction corrections, it should be preserved even in the presence 
of final-state interactions. Interestingly, the same high-energy behavior of the eigenphase 
shifts that makes the elastic scattering amplitudes purely imaginary leads to the almost real 
two-body decay amplitudes. 



II. BASIC PROPERTIES OF THE INELASTIC DECAY AMPLITUDES 

A. Analyticity 



We consider the weak decay 



H — » hadrons (1) 



where if is a heavy particle such as the D and B mesons. The final state is generally a 
multiparticle state of N hadrons. Going off the H mass shell, we call the H mass squared 
as the variable s and examine the analytic property of the invariant decay amplitude M(s) 
in the complex s-plane. 

The analytic property of the S-matrix elements was extensively studied decades ago . 
We obtained many rules of computation by examining the Feynman diagrams though a 
rigorous proof without referring to the diagrams was given only in the limited cases. To 
avoid inessential complications, we consider the case where H and the final hadrons are all 
spinless. We shall use the in- out formalism to simplify our notations. The invariant 



amplitude for the decay into the A^-hadron state / is defined by 

{r\p,)\H^\HiP)) = M^ip,), {Y.p. = P) (2) 

i 

where the one-particle states are normalized as (p|p') = (27r)^5(p — p'). M^{pi) is actually 
the function of all possible Lorentz invariants made of Pi and P. Aside from s{= P^), there 
are (A^ -|- 1)(A^ — 2)/2 independent invariants, which we may choose as 

S^J = {P^+P,f {^>J{J^N-l)), (3) 
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since there is one linear dependency relation, 

s = 2Y^ip,-p,) + J2ml (4) 

i>j i 

We may include short-distance strong interactions in Huj by using the QCD-corrected effec- 
tive Hamiltonian. Though we shall treat as a local operator, the analytic property to be 
discussed below does not depend on locality of H^. What is important is that transfers 
no energy- momentum. 

To study analyticity, we introduce the auxiliary (unphysical) amplitude, 

{r{p,)\H^\H{P)) = M'{p,). (5) 

We can obtain M~{pi) by imposing the incoming boundary condition on the hadrons in 
Eq.(^. Diagramatically, it amounts to flipping the sign of ie in all Feynman propagators. 
If is time-reversal invariant, it holds that 

{r\p,)\H^\H{P)) = {H{P')\HM{pT). (6) 

where p'^ and P' are obtained from pi and P by reversing the signs of their space components. 
Eq.(^ reads 

M+(p,) = M-(pO*, (7) 

where the asterisk denotes a complex conjugate. The reversal of the signs of the space 
components of pi does not change the Lorentz invariants, Sij and s. Since both M^{pi) and 
M~{p'j) are real below all thresholds and therefore coincide with each other there, Eq.(|^ 
means that they actually represent values of a single analytic function on two different 
Riemannn sheets. In terms of the Lorentz invariants, Eq.(|^) can be written as 

M(s + ie, Sij + ie) = M{s - ie, Sij - ie)* . (8) 

When is not T-invariant, it is convenient to work with each short-distance-corrected 
weak Hamiltonian individually after separating out a T- violating phase; = H^e^^^^^ . 
Then Eq.(^ is vahd for H^. We shall use the words " T- violation" and " CP- violation" as 
equivalent, assuming the CPT invariance. 

One problem about the multibody decay amplitude is that there are so many variables; 
three independent variables even for the three-body decay, one more than in the Mandelstam 
representation for two-body scattering. More a serious obstacle is that the real analyticity 
relation Eq.(|) holds only when we go across all cuts in s and Sij from +ie to —ie simul- 
taneously. We are able to write a dispersion relation in one of the variables keeping the 
others fixed, for instance, in variable s keeping Sij{{ij) ^ {N, N — 1)) fixed to Sij + ie. Then 
M(s — ie,Sij + ie), which is not simply related to the physical amplitude, enters the dis- 
persion integral. In the case of a three-body final state, for instance, such an amplitude is 
the complex conjugate of the unphysical amplitude M(s -|- ie, Sij — ie) in which the particle 
pairs (1,2) and (1,3) interact with the wrong sign phases —S12 and — 5i3, respectively, while 
the particle pair (2,3) interacts with the right sign phase 523- Only for the two-body decay 
amplitudes does Eq.(H) give the simple real analyticity relation: 
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M{s + ie) = M{s-iey, 



(9) 



so that the discontinuity in the complex s-plane is directly related to a physical process. For 
the decays of > 3, we can write only multivariable dispersion relations as an extension of 
the Mandelstam representation. We see little chance of extracting a useful information out 
of them. 



B. Dispersion relation 

The standard dispersion relation relates the real part of M{s) to the imaginary part. To 
relate the phase to the magnitude, we write the dispersion relation for logarithm of M{s). 
Notice that In M{s) has the same analytic property as M{s) except for cuts due to zeros of 
M(s). Choosing the contour of the Cauchy integral as usual, we can write the dispersion 
relation in the once-subtracted form: 

In M(s) - In M(0) = - T -J^^l-ds', (10) 

TT J so S'[S' — S) 

where sq is the lowest threshold of final states. The phase 6{s) should be normalized to zero at 
s = So in this representation to keep In M(s) finite at s = sq. When Eq.(|TO|) is exponentiated, 
it is the representation of Mushkelishvili which was first applied by Omnes to a study of 
the electromagnetic form factors. Hereafter we shall refer to this exponentiated dispersion 
relation as the Omnes-Mushkelishvili representation The application was limited to the 
energy region where the two-body scattering is elastic. Contrary to some misconception 0, 
however, the same dispersion relation can be derived for the two-body decay even in the 
presence of inelastic channels. We emphasize this point since otherwise the representation 
will be of no use in the heavy hadron decay. The key observation here is that once the phase 
along the cut is known, a real analytic function is unique. 

Actually there is one uncertainty that cannot be fixed by the analyticity. It arises from 
possible zeros of M(s). If M(s) has zeros at Si {i = 1, 2, 3, ■ ■ ■), they generate logarithmic 
singularities for InM(s) and contribute to the dispersion integral. When such zeros are 
included, the amplitude M(s) is expressed as 

M(si^P(.)e.,{t££^,,), („) 

where P{s) = M(0)nj(l — s/si). We need some physical argument to determine the poly- 
nomial P{s). In the days of the bootstrap theory we used to resort to a certain philosophy 
of determinism: there must not be a free parameter which we cannot control on. We do not 
think that we can argue along the same line in the context of QCD. In our case we avoid 
this problem as follows. 

Let us write Eq.([TT|) back in the form 

InM(s) = \nP(s) + Mm -( -^^ds' - ^-^ds') (12) 

^ ' ^ ' A2^oo TT \Jso S'-S Jso S' J ^ ^ 

and then define M{s; m*^) by 
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lnM(s;m*2) =lnP(s)+ lim -( -^^ds' - T ^rfs'V (13) 

^ ' ' ^ ' A2^oo 71 Vim-2 S'-S Jso S' J ^ ^ 

M{s; m*^) is the unphysical decay amplitude in which the final-state interactions of the 
energy range from the threshold to ^/s = m* has been removed. In terms of this amplitude, 
the physical decay matrix element M{s) is expressed as 



M{s) = M{s;m*^)expi- r ^f^ds']. (14) 

VtT Jso S' — S / 

In this form all strong-interaction corrections below ^/s = m* are explicitly factored out 
in the exponent. The correction factor includes both short- and long-distance corrections 
except for the interaction responsible for formation of hadrons. Concerning the hadron 
formation interaction, we encounter a fundamental issue in the final-state interaction theory. 
Both the final-state interactions and the formation of hadrons result from the same QCD 
force. Nonetheless, in order to formulate the final-state interaction theory, we must separate 
the hadron formation forces from the long-distance interactions between hadrons.^ The 
representation above provides one way to do so. While M{s; m*^) cannot be measured in 
experiment, it is the closest to what theorists have been calculating for decay matrix elements 
by various methods without final-state interactions. The appropriate choice of a value for 
m* in Eq.(|T^) is 

niH < m* < mw- (15) 

A remark is in order on the two-body final state. The dispersion relation in Eq.(|T^ holds 
for any two-body final state. In the elastic energy region, we normally choose an isospin or 
an SU(3) eigenstate for it. The reason is that such a state is an eigenchannel of the strong 
interaction S-matrix, and therefore that the phase is identified with the strong interaction 
phase by Watson's theorem |jll| . At high energies where two-body states couple to multibody 
final states, however, two-body final states are no longer eigenstates of the S-matrix, no 
matter which isospin states we may choose. Then it happens that the net phase 6{s) of 
the two-body decay amplitude depends not only on strong interactions but also on weak 
interactions, even when one computes amplitudes due to a single effective weak Hamiltonian. 
Though this was already pointed out @] in the past, it is worth emphasizing since it is the 
origin of all complications when rescattering is inelastic. We shall see the point more clearly 
in Section III. From this viewpoint, two-body isospin eigenstates in the inelastic region are 
just as bad as the J3 eigenstates of indefinite isospin in the elastic region. When inelastic 
channels are open, writing for instance the decay amplitudes for D~tt'^,^ D 7r° in 

the J = 1/2 and 3/2 amplitudes contributes to very little to solving the problem. 



^Recall that the hadrons were considered as elementary when the phase theorem was proved |11| 



6 



C. High-energy behavior of the phase 



Imagine that 6{s) approaches the asymptotic value 6{oo) at some value m well below 

6{s) ^ 6{oo) {m < y/s < mw)- (16) 
The exponent of Eq.(|l^) can be estimated for such 5{s) as 

1 Ks') ^ 1 ( r- ^ r\ 

TT J so S' — mjj — le IT J so / s' — mjj 

= ^ ' In +Z7r + ^ . 

The contribution of 0(m^/m|^) from the region below ^/s = m is negligibly small if 
m'jj 3> m^. The dominant contribution comes from the asymptotic energy region. When ex- 
ponentiated, this integral generates an enhancement or suppression factor of (m* /mHy^^°°^^^ 
for M{s). If, for instance, the phase reaches 6{oo) = ±7r/2, the strong interaction correction 
would alter the amplitude by factor 16 for rriH = ttlb and m* = ^m^r. It means an enhance- 
ment or a suppression of factor ~ 250 in rate. There is no evidence for such a huge dynamical 
enhancement or suppression when we compare the observed two-body decay rates with the 
theoretical estimates in the D and B decays. The so-called color suppression observed in the 
B decay should be attributed not to a severe dynamical suppression by strong interactions 
but to lack of strong-interaction corrections, since the suppression exists prior to final-state 
interactions. 

The obvious alternative to the asymptotic behavior of Eq. ( ]T6| ) is 

6{s) im<^/s< mw). (18) 

The enhancement or suppression is milder in this case. 

The fractional power {m* /mHY^^°°^^'^ does not appear in the conventional calculation of 
the decay matrix elements. The short-distance QCD corrections enter in fractional powers 
of (lnm*/lnm/^). //we wish to reproduce the short-distance correction factor of the renor- 
malization group with our final-state interaction integral, we would choose such that 6{s) 
approach zero asymptotically as 

6{s) — > 6/ Ins, (m < < mw) (19) 

where 5 = 7i„/6o, 7«) is a constant determined by the anomalous dimension of Hy^, and 
bo is from the running QCD coupling. The asymptotic behavior of Eq.(|19]) leads to 



^Since there is no large characteristic energy scale of strong interactions, it is reasonable to assume 
that M(s) approaches its asymptopia as early as hadron scattering amplitudes do. If 6{s) should 
keep oscillating, InM(s) would behave like e''*' along some direction in the complex s-plane, which 
would prevent us from writing the dispersion relation to start with. 
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(Inm* / InmHy^^^ for the amplitude. Since our correction factor includes both short- and 
long-distance effects, Eg . ([T9|) should represent only the portion of the asymptotic phase that 
is attributed to the short- distance interactions. We must study the long-distance effects by 
some other means. 

The main result of our dispersion relation is that the phase of the two-body decay 
amplitude of a heavy hadron should be zero or a small value for a large initial mass {s —>■ oo). 
It should certainly not be ±90°. It is worth noting that, unlike the phase, the magnitude of 
the amplitude can be subject to substantial final-state interaction corrections since it picks 
up the effects of all energies. 



III. EIGENPHASE SHIFTS 

In this Section we study the decay phases from the viewpoint of the eigenphases of 
the S-matrix of high-energy scattering. We need to know about the composition of the 
eigenchannels and their high-energy asymptotic behavior. We shall use the dual resonance 
model as our guide since it is the most successful model that incorporates the relevant 
aspects of the quark model spectroscopy and long-distance hadron scattering. 

The partial- wave S-matrix of strong-interaction is diagonalized in terms of the eigenchan- 
nels carrying the quantum numbers of hadron H. Labeling the eigenchannels by |a), ■ • • 
we can express the S-matrix elements as 

Sab = (&°"*|a^") (20) 

Without loss of generality we normalize the eigenphase shifts to zero at their respective 
thresholds. Experimentally observed are the hadronic states with each particle carrying 
definite momentum. We project those states onto the eigenstates and denote them by 
h. By completeness of eigenchannels, h can be expanded as 

\h} = J2o,M, (21) 

a 

where we can choose Oha to be an orthogonal matrix as a consequence of T-invariance for 
strong interactions. 

The partial- wave amplitude a"^(s) for elastic scattering is expanded in the eigenphase 
shifts as 

a'{s) = Y.0ly'^sm6a, (22) 

a 

or in the real and imaginary parts, 

Rea'^(s) = Ol^ cos 6a sin 6a, (23) 

a 

lma-^{s) = J20l,sm^6a. 

a 

Both Oha and 6a are s-dependent. Strictly speaking, once multibody states are included, we 
must label the eigenchannels with continuous parameters. Therefore the discrete summation 
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in Eqs.(pTD and ( p4D is symbolic. The state density per unit energy in a volume characteristic 
of strong interactions (~ m~^) may substitute as the effective number of states for multibody 
channels. 



Actually the dual resonance model ||T2[ of the late 1960 's answers to how many states 



exist at energy ^/s. In this model, the number of states was counted to be [jl3 



^0 ~ (^.^)(rf+i)/2 exp(v^/mo), (24) 

where a' is the Regge slope (~lGeV~^), d is the space-time dimension, and mo = 
{3/2a'dY^'^ /n, which takes a value ~ 200MeV for d = 4. This state density uq contains the 
states of all angular momentum J < (a'sY^"^ at mass ^/s for the interval of ~ lGeV~^(~ a'). 
While the highest J state has no degeneracy, degeneracy of states rapidly increases with de- 
scending J. It also includes the states of negative norm on the daughter Regge trajectories 



in the case of 4-dimensional space-time. Nonetheless the state density of Eg. (|24D , particu- 



larly the exponential dependence, gives an order-of-magnitude estimate for the number of 



the J=0 states in the dual resonance model. Hagedorn [|1J] introduced a statistical model of 
hadrons with quite a different motivation. It is amusing that his model led to essentially the 
same state density with a very close value (~160 MeV) for ttlq but with a slightly difference 
power of s in front. In Hagedorn's model, the closeness of mo to the pion mass was explained 
by the fact that every time energy increases by mo, one more pion evaporates and causes 
the exponential growth of the state density. 

The phenomenological success of the dual resonance model confirmed that the s- 
channel resonances are dual to the non-Pomeron Regge exchanges while the Pomeron is 
dual to the nonresonant continuum in the s-channel The dual resonance model at tree- 
level incorporates only the non-Pomeron trajectories. The Pomeron term was included by 
adding the two-body nonresonant intermediate states. At high energies, the resonant states 
are so broad in width that they are not recognized as resonances but make up the smooth 
Regge asymptotic behavior of non-Pomeron exchange. In the yet higher energy region, 
the diffractive scattering dominates so that the s-channel states consist almost entirely of 
nonresonant states. Henceforth we shall call the s-channel states dual to the Pomeron and to 
the non-Pomerons as the nonresonant and the resonant channels, respectively, even though 
no resonance peak appears in the resonant channels at high energies. We shall also use the 
words, the diffractive and nondiffr active channels, for them. In the dual resonance model, 
the effective number of the nonresonant channels is even larger than tiq of the resonant 
channels in the high-energy limit. 

The distinction between the Pomeron and the non-Pomeron trajectories is best described 
by the quark diagram (See Figure 1.) The non-Pomeron exchange in the boson- 



boson scattering is described by a pair of quark-antiquark in the intermediate state (Fig. la) 
which represents a tower of resonances in the s-channel and at the same time the Regge 
trajectory exchanges in the t-channel. In contrast, the Pomeron exchange corresponds to 
the "disconnected" quark diagram of the qqqq four-quark intermediate states in the s-channel 
as shown in Fig. lb. In the context of QCD one pair of qq exchanges relatively soft gluons 



with the other pair of qq in the Pomeron exchange |18|. This quark diagram not only ensures 



the exchange degeneracy for a pair of non-Pomeron trajectories with opposite signatures. 



as observed in scattering experiment |T^, but also explains the absence of the mesons with 



exotic quantum numbers such as /> 3/2 in meson spectroscopy. 
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A. Final-state interactions in eigenphase shifts 



The T-invariance relation in Section II can be easily extended to the decay amplitude 
into the eigenchannels. Using the completeness condition J2a = 1 and (a^^'la*") = 

g2«(5a^ obtain 

{a"^'\H^\H) = e^'^^^'\a'"''\H^\H)\ (25) 

where s = mjj. Accordingly the decay amplitude Ma{s) into the eigenchannel a carries the 
eigenphase Sa{s). Inserting the complete set of eigenchannels and using the expansion of 
\h), we obtain 

{h"^'\H^\H) = ^Of,a{a"'''\H^\H). (26) 

a 

The decay amplitude Mh{s) for H h is now expressed in the terms of the eigenchannel 
decay amplitudes Ma{s): 

Mh{s) = Y.OhaMa{s). (27) 

a 

Separating the phase 6a{s) from Ma{s) as 

Ma{s) = Ma{s)e''^'^'\ (28) 

we put Eq.([27|) in the form 



Mh{s) = Y.OhaMa{s)e''^^'\ (29) 

a 

Comparing Eq.(p9D with Eg. (p2|) , we clearly see that the net phase of Mh{s) has little to do 
with that of a'^(s). The phase of Mh{s) agrees with that of a'^(s), barring an accident, only 
when there is only one eigenchannel so that the elastic unitarity holds, 

11 + 2^-^(5)1 = 1. (30) 

It was pointed out that the partial-wave projection of the diffraction amplitude for the tttt 
scattering at the B mass is far short of the unitarity limit 0]; |1 + 2m'^(s)| ~ 0.6. For the 
D decay, which occurs near the resonance region or a little above it, the nonleading Regge 
terms cannot be ignored in two-body scattering. We shall see in Section V that even after 
adding the non-Pomeron terms the elastic unitarity does not hold at the D mass. 

The amplitude Ma{s) defined in Eq. (|28|) still contains strong-interaction effects. To 
define the eigenchannel decay amplitude free of final-state interactions, we must separate 
out not only the phase correction but also the magnitude correction Aa(s) by Ma{s) = 
Moa{s)Aa{s). If all inelastic channels are approximated as two-body or quasi-two-body 
states, Aa(s) can be written in the Omnes-Mushkelishvili representation: 

A.W=expf^r''^^.,'l (31) 



71 Jsoa S' — S 
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where V stands for the principal value integral. One obvious property of Aa(s) is that it 
is positive definite. Whether Aa(m|^) gives an enhancement (> 1) or a suppression (< 1) 
depends on the sign of 5a{s) for all values of s from the threshold Soa to m^, not just the 
on-shell value 5{m'\j). In the nonresonant channels, we shall argue later that 6a is small in 
magnitude and the sign of 5a can be easily flipped as energy changes when a large number 
of channels mix. If so, a correlation of the magnitude correction Aa{rn\j) with the on-shell 
phase value Saimfj) is tenuous, if any. With the magnitude enhancement factor written out, 
Eq.(^9]) turns into 

M,(s) = Y.OhaMoa{s)Aa{s)e''^^'\ (32) 

a 

For the resonant channels, one can compute the Feynman diagram for the decay process 
H i?(resonance) /i(two-body) through a resonance R. The decay amplitude for the 
resonant eigenchannel r takes the form. 



Mr{s) = fHR , — , (33) 

m% - s - imRltot 

where fuR is the H-R pole transition strength and Ttot,r are the total and partial decay 
widths of R. Note that Vtot = O(mj^). In the quark model Jhr is the overlap of the wave 
functions of H and R. Since the overlap does not increase with s, M^(s) decreases like 1/s 
or faster as s ^ oo. If we express the resonant channel contributions in the form of Eq.(^), 
it means Mor(s)Ar(s) ^ 1/s or faster for each resonant channel r. 



B. Strong interaction scattering 

Elastic scattering provides useful pieces of information about the eigenphase shifts. Ex- 
periment shows that the imaginary part dominates over the real part for the invariant 
amplitude T{s,t) of elastic scattering at high energies. Theoretically, the dominance of 
lmT{s,t) is a general consequence of analyticity and crossing symmetry, not specific to the 
Regge theory, when the total cross section approaches a constant up to powers of logs. The 
amplitude of the fiat Pomeron trajectory, 

T{s,t) = tsatote'', (34) 

gives a reasonably good description of the diffractive scattering in the whole high-energy 
region relevant to us.f\ Eq. (p4D leads to the elastic cross section,^ 



^At energies above ^/s ~ 100 GeV, the total cross sections actually rise very slowly with energy. 
One fit to pp-collisions gives atot{s) ~ I3p{s / sof''^^ + /3p_/(s/so)~'^'^^ This s-dependence 

requires that the forward scattering amplitude contains a real part: Rer(s, 0)/ImT(s, 0) 
tan(7re/2), where e = 0.08. The forward scattering amplitude contains a real part by about 
10% even at extremely high energies. However such high energies have no direct relevance to the 
final-state interaction of the B decay. 

^We ignore all hadron masses as compared with ^/s throughout this section. 



11 



(^ei = crlJ327ib, 



(35) 



and to the partial-wave amplitude 

a-'is) =i(Ttot/l6Trb, (J < (s/so)'/'). (36) 

When we parametrize T{s,t) by Eg. (p4D , the partial-wave elasticity crii/cyiot (= Ima'^(s)) 
for J <^ {s/sqY/'^ is equal to twice the total elasticity, {1/atot) J{daei/dt)dt. The values for 
(Jtot and b can be extracted from the experimental data on pp, np, and Kp scattering. The 
factorization of the Regge residues]^ relates the high-energy total cross sections by 

^MM' _ Mp M'p I pp ,or7\ 
'-'tot — tot '-'tot I'^toti I -J 'J 

where M and M' stand for mesons. With crfo^ = 37mb, cr^j^ = 21mb, and a^f = 17mb for 
the diffractive contribution of atot at ^/s = 2 ~ 8 GeV [|19|, we obtain 

<l = 12mb a,^,- = 10mb. (38) 

The numerical values are roughly in line with the empirical law of the quark number counting 
for the total cross sections: (jfot'^' : (jf^t '■ '^lot = 2^ : 2 x 3 : 3^. The diffraction width 
parameter b obeys the inequahty bpp > b-^p > bxp- In one analysis [ pO(] 

6pp~5GeV~^ 6^p~4.3GeV"^ ~ 3.2GeV-l (39) 

The parameter b is related to the effective target size of colliding hadrons in elastic scattering. 
The inequality bpp > b-j^p indicates that the proton is a little more spread than the pion. The 
electromagnetic form factors show the same trend: (?"p/6)^/^ ~ ^J2m~^ and (r^/6)^/^ ~ ^• 
The relation bj^p < b^^p may be interpreted as a result of the less intense interaction of the 
s-quark with the u/d-quarks, which reduces the effective size of K as well as afof. This line 
of argument leads us to the 6-parameters for tttt and Kti scattering somewhere around 

b^^ ~ 3.7GeV-^ bKn ~ 2.8GeV"l (40) 

If we extrapolate this reasoning to the Dtt scattering and ignore the c-quark interaction with 
the u/d-quarks, we are led with a^J = and fo^Tr < bxn to 

af^J = 6mb bon ^ 2.4GeV~l (41) 

When we substitute these values of parameters in the partial- wave projection of Eq.(^6l), we 
obtain 

r 0.16 (tttt) 

Ima^(s) - < 0.17 (Kvr) (42) 
[ 0.12 (Dtt) 



^ The factorization can be proved only for relatively simple J-plane singular lities. It is an assump- 
tion for the more general Pomeron. 
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The value of Ima'^(s) is even smaller for ipn and ipK. The precise values of the right-hand 
side of Eq. (^) are not important in the following. It is no surprise that the elastic unitarity 
is not satisfied for any of the above processes: 

|1 + 2m^(s)| = 0.66 ~ 0.76. (43) 

Note that |1 + 2ia'^{s)\ = 0.5 is the limit of the completely absorptive black target. When 
a"^(s) is purely imaginary, the value of |1 + 2ia'^{s)\ does not tell the whole story. The 
partial- wave inelasticity reveals more of the scattering mechanism. The numbers in Eq. (^2p 
give the partial-wave inelasticity (=1— Ima^^) in the range of 

^Lei/^L = 0.83 ~ 0.88. (44) 

It should be pointed out here that the s-wave phase shift is dominated by long-distance 
physics. The largest contribution to a'^(s) for small J comes from the region of the mo- 
mentum transfer —l/b<t<0. The contribution from the perturbative QCD region of 
large \t\ is small. It is ~ l/|t|^/^ that determines the effective distance of interactions. It is 
true, however, that the s-wave contains a larger share of short-distance physics than high 
partial-waves. 

Comparing Eq. (|4^) with Ima"^ of Eq. (^21) , we find that the average or typical eigenphase 
shift of the diffractive channels should be in the range of 

sin^^rf ~ 1/8 ~ 1/6, (45) 

where the subscript d for the phase shift stands for "diffractive" . The smallness of Rea'^/Ima'^ 
for the diffractive scattering suggests that 5^ (c? = 1, 2 ■ ■ ■) spreads over positive and negative 
values (modulo nvr) in an approximately symmetric distribution with respect to 5^ = 0. For 
such a distribution a large cancellation occurs among different eigenchannels in Rea'^(cx 
cos(5rfsin5rf) while every term adds up in Ima"^(oc sin^ 6d)- What about the rni ambiguity for 
5d? Since the diffractive channels are nonresonant, 6d starts at and does not go over 7r/2. 



Because of Wigner's theorem |2^ on the causality constraint on phase shifts, it is not very 
likely for 6d to turn clockwise and go over —n/2. Therefore 6dS stay in —tt/2 < 6d < 7r/2, 
spreading symmetrically with respect to 5^ = over the range roughly 

- 7r/8 <Sd< n/8. (46) 

Here again, the precise values of the upper and lower bounds are unimportant to us. 

Even in the decay into a diffractive multibody channel a resonance can appear in a 
subchannel giving rise to a large phase when the subchannel invariant mass coincides with 
the resonance mass. For instance, H — > i?7r (nonresonant )—*> KKn at rrij^j^ = rriji. If R 
is a sharp resonance of KK, we treat the process as a two-body decay into Rtc. If not, 
the final-state is a three-body channel. Though the decay amplitude has a large phase at 
rrij^-j^ = niR, this phase is washed out after the rrij^j^ is integrated over with -y/i fixed to niH- 

Let us turn to the nondiffr active channels. Their contribution to atot falls off like l/s^^^ 
or faster relative to that of the diffractive channels. To obtain the subdominant terms in 
s, we add the nonleading Regge contributions. One may wonder about a possible isolated 
s-channel non-Regge singularity which may contribute only to a single angular momentum. 
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Especially relevant is the fixed singularity a'^ oc Sjo/ J. This singularity generates a constant 
term in T{s,t) for all values of s and t. Since hadrons are composites of the scale Aqcd, 
there must not be such a hard interaction. 

According to the Regge duality, the nondiffr active and diffractive portions of the am- 
plitude are separately dual to the resonant and nonresonant states, respectively, in the 
s-channel. Therefore the nondiffr active or resonant contribution a;!(s) to the partial- wave 
amplitude can be expressed as 

a'As) =^0W^- sin 6r, (47) 

r 

where the subscript r stands for "resonant". In these channels, the eigenphase shift 5r turns 
around counter clockwise slowly passing 7r/2 at the resonance and approaches vr asymptoti- 
cally. Comparison of the asymptotic energy dependences of the diffractive and nondiffractive 
amplitudes gives us 

where Ppp^f are the properly normalized Regge residues at t = of the Pomeron and the 
P-/2 trajectory, and s is in the unit of 1 GeV^. We can extract the Regge residues from the 



energy dependence of total cross sections. For n^p (vr p) scattering for instance. 



/3^_^//?p^ 0.75(1.18). (49) 

We can relate the left-hand side of Eq. ( ^HD to the numbers of the resonant and nonresonant 
eigenchannels. By replacing sin^ Sr by unity for the resonant channels and substituting the 
average value for sin^ 6d, we obtain 



X^r sin^ 5r(s) ^ rir 



(50) 



With Eqs.(E|) and (M), we obatin 



nr^n{sm^6d)/s^^\ (51) 



which agrees with our expectation rij. <^ n. 



IV. RANDOM S-MATRIX APPROXIMATION 

When a very large number of eigenchannels are present, studying individual channels is 
impractical. To study the spacings and widths for hundreds of the densely populated reso- 
nances in complex nuclei, nuclear physicists introduced a certain randomness hypothesis in 
the multichannel S-matrix. The work was started by Wigner |23|, pursued by many p5[, and 



brought into a mathematical sophistication by Dyson |26|. It succeeded in reproducing vari- 
ous features of those resonances |^ . A similarity of the nuclear resonances to the multitude 
of the hadron channels in our problem suggests us to study physics of the uncontrollably 
many eigenchannels with the randomness hypothesis. 



14 



The randomness of the channel mixing Oha can arise in our problem when a very large 
number of eigenchannels exist in degeneracy with the two-body final-states. Since even 
a small coupling strongly mixes a pair of degenerate states, any two-body state in the 
highly inelastic energy region is a linear combination of many eigenchannels. The expansion 
coefficients of a given state into eigenchannels, namely Oha, is sensitive to the strength of 
channel couplings, but many features of physics, for instance 6a, should not be sensitive to 
small variations of channel couplings. It is therefore reasonable to postulate that quantities 
of our interest can be computed by replacing products of Oha with their statistical averages 
over the phase space of the 0{n) rotations of Oha- After the 0{n) average has been taken 
for products of Oha, we are left with 6a- In our problem unlike the nuclear resonances, we 
cannot postulate a complete randomness for the distribution of 6a- If we assumed that both 
6a and Oha are completely random or chaotic, Eq.(plD would lead us to {OhaOkb) = 6ab6hk/n, 
(sin^ (5a) = 1/2, and (cos(5aSin(5a) = 0. Consequently 

Rea-^ = 0, Ima-^ = 0.5. (52) 

This corresponds to the scattering from a black disc that gives the 50% elasticity due to 
the shadow scattering in disagreement with Eq. (P2|) . In the actual high-energy scattering, 
o'ei/o'tot is considerably less than 50%. It means that a hadron target behaves like an opaque 
disc at high energies. To describe such a target, we must postulate randomness only within 
the restricted range as specified in Eq.(^6D.p| 



A. Dominance of the real part for decay amplitudes 

We start with Eq.(P2|) and isolate all eigenchannel dependences by expressing Ma{s) = 
Moa^a^^^" back in terms of the decay amplitudes of the hadron basis k which are free of the 
final-state interaction (of < mw)- We denote such decay amplitudes by Mofe(s). Writing 
the nonresonant and resonant channels separately, we have 

Mokis) = OkdModis) + E OkrMoris). (53) 

d r 

Substituting the inverted relation of this in Eq. (|52D , we obtain 

Mh{s) = ^0,,Ofc,Mofc(s)Arf(s)e*'^'*(^) + ^0;,,Ofc,Mofc(s)A,(s)e^^'-W. (54) 

d,k r,k 

Using {OhaOkb) = 6hk6ab/n, we obtain the random phase values for the real and imaginary 
parts of M/i(s): 

Mh{s) = Moh{s){{A cos 6) + i{A sin 6)), (55) 



^We previously studied mainly the magnitude of the squared decay amplitude in this statistical 
model Here we focus on the decay phase by refining some of the previous postulates. In Ref 
0, the magnitude factor Aa(s) was not separated but ignored for simplicity. 
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where 



(Acos5) = - (Y.^d{s) cos 6d{s) + Y.Ar{s) cos 6r{s)), (56) 

and cos 5 sin 5 for (A sin 5). 

Let us leave out the resonant eigenchannels for the moment. Then, first of all, the phase 
of Mh{s) approaches a common limit in Eq.(^) independent of isospins or charge states of h, 
say, K'Tc'^ and TC^tt^. The common limit does not depend on the effective weak interaction 
H^. This marked simplicity is valid only in the random phase limit. Considering that Ima'^ 
takes roughly the same value for all meson-meson scatterings (cf. Eg. (^2]) ) , we expect the 
common decay phase is not very sensitive to the final hadrons. 

We have pointed out earlier that in the quasi-two-body approximation to the inelastic 
channels, Aa{s) is positive definite and uncorrelated with the on-shell Sa{s). If this is the 
case, the terms of different eigenchannels cancel each other in lmMh{s) because of the 
random signs of sin 6ci:f\ 

ImMhi-s) = Ad{s) sinSd - E Ad'{s)\smSd'\. (57) 

Sd>0 <5d/<0 

In contrast, all eigenchannels add up in the real part since cos 5d > in the restricted range 
of Eg. d^) . It is similar to the situation in the elastic scattering amplitude of Eg. , but 
this time the cancellation occurs in the imaginary part instead of the real part. Because 
of the deviation of Arf(s) from unity due to enhancement and suppression, the cancellation 
may not be as good as in the scattering. We can set with Egs. (p5D and ( p6D a loose upper 
bound on the imaginary-to-real ratio for Mh{s). In terms of tanSh = lmMh{s) /R.eMh{s) , it 
is given by^ 

I tan^ftl < |(sin(5)|/(cos5). (58) 

The right-hand side is less than ~0.4 acccording to Eg.(^). This number would be realized 
when no cancellation occurs in Eg. {^7\) . Since Ad{s) and the sign of the on-shell 6d{s) are 
only tenuously correlated, we expect in reality a fairly high degree of cancellation between 
the terms of 5^ > and 6d < 0. Therefore the actual value of the imaginary-to-real ratio is 
most likely much smaller than 0.4. By pushing our approximation further, let us set A^ to 
a common number and cos^d ~ 1. Then we obtain 

tan(5/i ~ ^^^(sin^ (5). (59) 
Ima-' 



^The possibility of many phases averaging out to a small net decay phase was mentioned earlier 
by Wolfenstein to the author. 

Section II, 5h was simply written as 6{s). 

^If the sign of the on-shell 6d is not correlated with A^^ > 1 or < 1, setting A^^ to a number 
independent is justifiable. 
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The right-hand side is zero for the purely diffractive amphtude. Even if we use the 
parametrization atot ~ s*^'^^ fitted to much higher energies which requires a nonneghgible 
real part for T(s,0), the right-hand side is ~ 0.047r x {siv? 5) ~ 2 x 10~^. 

Next we look into the contribution of the resonant channels to the decay phase. It 
falls like s"''"^"^ s^^/^ relative to the diffractive contribution. In addition to the s"^/^ 
suppression, the chiral structure of the weak interaction suppresses the transition from the 
initial state H oi = 0~ to the resonant channels. Let us examine this suppression. 

Keeping in mind that in the dual resonance model the resonant channels are made of 
two-quark states in the quark diagram (Fig. la), we examine the transition of H{= qQ) to 
the resonant state R{= Qq')^ The bare weak Hamiltonian which causes the transition is 
the four-quark operator, 

= 4(G/v^)l^;,r,Q(7^g^)(g^Qi), (60) 

where we have suppressed the Dirac structure of the quark bilinears. The short- distance 
QCD corrections to if^; induce the effective interactions such as (5''L(Aa/2)gj^)(g^(Aa/2)(5i) 
and the penguin interaction {q^{\a/2)qii){q'i^{Xa/2)QL). The weak transition H ^ R is the 
so-called annihilation or exchange process. The matrix elements of the annihilation and 
exchange processes are suppressed by the decay constant ///(= O^f.,,)) of the meson H and, 
for the non-penguin interactions, also by the factor (m^/ + niq) of chirality mismatch on the 
side of R. The suppression due to the first one, fn/mH, together with the energy dependence 
suppression s~^/^ or n^/n is severe enough to make the phase difference contribution of 
the resonant channels negligibly small. It should be reminded that the resonances in the 
resonant channels are so broad {Th = 0{mii)) that no sharp resonance enhancement arises. 
The suppression of fn / itlh occurs for the transition H R through the penguin interaction 
too. 

In short, the dominant decay process is through the spectator diagrams which lead to 
qqqq not qq in the final state. The decay H qqqq followed by a pair annihilation into qq is 
no other than the annihilation or exchange process, as we can see by drawing the diagram. 
The rescattering of qqqq without a qq annihilation is a diffactive process not a resonant one 
in the sense of the dual resonance model. Therefore, we can conclude that the decay phase 
due to the resonant channels are negligibly small. We estimate a typical magnitude of the 
phase differences of this origin as 

A5^ifH/mH){lGeV/s'/^) (61) 

for the decay modes where the spectator process dominates. When the chiral mismatch 
occurs for R, a decay matrix element is suppressed by the additional factor of iriqi/mH- 

B. Phase difference between amplitudes 

What we can measure in experiment is not the absolute phases but the phase differ- 
ences. Though the phase of {Mh{s)) is independent of the charge or isospin states of h, 



^'^The following argument is valid also for H — > R{= q"q') with minor modifications. 
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the fluctuation around the average value can be isospin dependent and therefore generate 
phase differences. The physical origin of this type of phase differences can be explained as 
follows: A pair of decay amplitudes have different compositions of diffractive eigenchannels. 
Therefore the interference between eigenchannels sums up to different net phases for the 
amplitudes. The phases of this origin are washed out in the random limit. The interference 
between the diffractive and nondiffractive eigenchannels is another possible source of the 
relative phases. 

We first examine the relative phase due to the fluctuations. The fluctuation is the 
standard deviation from the random phase limit of {Mh{s)). Averaging out the product of 
four Oha in the phase space of 0(n) group|^, we obtain for n ^ rir 

\AlmMh{s)\^ = {{ImMhis) - {lmMh{s)))^) , (62) 

= {\Mo\'){{A'sm'6) - (Asin(5))2) + 0(^|M,(s)|2), 
|AReMfe(s)p = (sin5 ^ cos 5), 

where the brackets denote the averages such as 

{\Mo\') = -T.\Mod\'. (63) 

d 

The sum J2a={d,r) \Moa\'^/2mH is the total rate for the decays induced by without strong 
interaction corrections. Since the decay amplitude is dominantly real, the phase difference 
is more sensitive to AlmMh{s) than to AReM/i(s). The relative magnitude of the standard 
deviation AMh{s) to the average |M/i(s)| depends on the channel number n. If n is nearly 
as large as the dual resonance model indicates, or even a small fraction of it, |AM/i(s)| 
is negligibly small and no relative phase arises from fluctuations. In this case, the two- 
body channel couples to an enormous number of inelastic nonresonant channels though its 
couplings to individual channels are very weak accordingly. The fluctuations almost even 
out when this happens. On the other hand, if n is so small that 1/n is comparable to the 
branching fraction to the channel h, AlmMft(s) can be close to |ImM/i(s)| itself. When 
|M/i(s)| is small by an accidental cancellation of high degree among different eigenchannels, 
the fluctuation can be large in proportion. To ensure that the random phase is a good 
approximation, therefore, we should apply its predictions only to the two-body modes having 
relatively large branching fractions. 

The phase differences arising from an interference with the resonant eigenchannels are 
suppressed by the chiral structure of the weak interaction operators and decrease like 
with energy. Their contribution is of the order of A6 as given in Eq. (|6l|) in general. How- 
ever, our argument fails in the modes where the spectator process is forbidden or highly 
suppressed. In such processes the diffractive eigenchannels dual to the Pomeron are missing 
or highly suppressed so that large phases and phase differences may potentially arise. It is 
unfortunate that the decay branching is generally very small for them. 



^^It is straightforward to obtain {OkaOiaOmhOnb) = {5ki5mn + 5km5in + ^kn^im) /■n{n + 2) for a = 6 
and = [{n + l)5ki5mn - hn^lm - 5kn8lm]l{n - l)n(n + 2) for a / 6 [§. 
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V. THE D AND B DECAYS 



We shall look into the specific cases of the D and B mesons with a few critical remarks on 
some of the recent attempts to compute the final-state interaction phases of the two-body 
decays of D and B. 

A. The D meson 

The most prominent two-body decay of the D meson is the Kir modes. From the observed 
decay rates of K 7i~^ and — > 7i^K~, K 7r°, the phase and magnitude of the decay 

amplitude were determined for D ttK of definite isospin. They discovered a large relative 
phase between the 1=1/2 and 3/2 decay amplitudes 

S3/2{ml) - 6^/2{ml) = (96 ± 13)° (64) 

with the amplitude ratio 

IM3/2/M1/2I = 0.27 ± 0.03. (65) 

Since the D meson mass is too low for the Pomeron to dominate in the Kn scattering, our 
argument in the preceding Section does not apply with a good accuracy. The phases of the 
s-wave Kn elastic scattering amplitudes ai{s) at energy niD was computed with the Regge 
theory and identified with the phases of the decay amplitudes for the Kn isospin eigenstates. 

The result of calculation in Ref. happens to be in a good agreement with the observed 
phase difference, Eq. (|6^) , within the uncertainties in the values of the Regge parameters. We 
argue however that such an agreement is fortuitous. In order to identify the scattering phases 
with the decay phases, the scattering must be elastic. The authors of Ref ^ computed only 
the phase difference of the Kn amplitudes, not the magnitudes of them. If they had done 
so, they would have obtained 

ai/2(m^) = 0.08 + 0.21i (66) 
03/2(771^) = -0.35 + 0.18Z. 

They do not satisfy the elastic unitarity: |1 + 2zai/2| = 0.60 and |1 + 2za3/2| = 0.95. 
The J = 1/2 amplitude is deep inside the Argand diagram while the 1 = 3/2 amplitude 
nearly satisfies the elastic unitarity. Actually, if one insists that Kn scattering be elastic at 
^/s = rriu, one can make 03/2 unitary in a model-independent way. Since the 7 = 3/2 channel 
is exotic, i.e., contains no s-channel resonances dual to the non-Pomeron trajectories, the 
imaginary part of 03/2 in Eq.(^) comes entirely from the Pomeron. The non-Pomeron 
trajectories can contribute only to the real part in the case of 03/2. Therefore one should 
make up for the unitarity violation of 03/2 by adjusting its real part: 

a3/2(m^) = -0.38 + 0.18z, (67) 

which gives 

6^/2iml) = 155°. (68) 
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Unlike 03/2 the unitarity deficit of ai/2 is very large. To recover the elastic unitarity for 
ai/2, the Pomeron daughters and the exchange-degenerate pairs of trajectories such as the 
daughters of P-/2 and of K*-K2 must contribute just as much as or more than the parent 
trajectories. If so, the Regge expansion itself would be unjustified at this energy. 

In the J = 1/2 channel there is an obvious candidate for a cause of inelasticity, that is, 
Kt] channel which should enter according to the flavor SU(3) symmetry. Then the {Kt^)i=i/2 
involves at least two eigenchannels, the 7 = 1/2 channels of the 8 and 27 representations of 
SU(3). In addition, there are multibody channels. The D meson decays into Kmrir of both 
1 = 1/2 and 3/2 including the modes such as K p and Kai. The combined branching frac- 
tion of all Ktctctt modes is even larger than that of Kit. Many = O"*" states are contained 
in KirTTTT, with which Kir mixes in general. It is not surprising if many eigenchannels are 
open for both 7 = 1/2 and 3/2 already at the energy of the D meson 

Is there any chance to identify the D —>■ Kir decay phase with the Ktt scattering phase 
all the way up to high energies despite the presence of inelastic channels ? We do not think 
so. In Eq. (|67D the real part of aj is positive for 7=1/2 and negative for 7 = 3/2 at the D 
meson mass. If it happens that ai/2 and 03/2 approach purely imaginary numbers at high 
energies, 61/2(3) would turn counter clockwise to tt/2 and 63/2(3) would turn clockwise to 
— 7r/2. The phase-amplitude dispersion of Section II then requires that Mi/2{D — ^ Kir) is 
enhanced over M3/2{D —>■ Kn) by factor [m* /mu)"^, where m* is the highest energy up to 
which one is willing to make the identification of the decay phase with the scattering phase. 
This relative enhancement factor is far too large unless m* is chosen to be comparable to 
rriD- 

To summarize, the D- meson mass is not low enough to allow identification of the final- 
interaction phase with the elastic scattering phase. Yet it is not high enough for our random 
S-matrix approximation to work with a good accuracy. The large phase difference obtained 
by experimentalists can be accommodated but not predicted. 



B. The B decay 

The meson-meson scattering at energy ^/s = is most likely asymptotic for vrvr, Kn, 
Dn or even for DDg. The resonances exist only at much lower energies, which means that 
the contribution of the nondiffr active scattering, i.e., the non-Pomeron exchanges, to the 
asymptotic scattering amplitude is negligible. For instance, if we make an estimate for the 
elastic Kn scattering at energy using the Regge parametrization of Ref [§, the diffractive 
contribution gives 

|l + 2m/| = 0.64 (7 = 1/2,3/2). (69) 

The nondiffractive contributions modify it only slightly into |1 + 2mi/2| = 0.63 and |1 + 
2203/2! = 0.65. The two-body scattering is clearly in the asymptotic energy region. The 
s-wave inelasticity for Kn is 

(Jinel/cTtot = 0.83 (70) 

from our estimate in Eq. (|42D . There is no chance for the elastic approximation to work at 
this energy. The fact that the largest branching fractions observed in the B decay are at 
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the level of 1% is a clear evidence for the presence of very many open channels. The B 
decay is the process where our results of Section IV should apply best. We shall restate our 
predictions for the B decay below. 

First of all, the two-body B decay amplitudes should be dominantly real up to CP viola- 
tion phases. According to our very crude estimate made in Eg. (pQ]) , the decay phases should 
be 7° or smaller. The relative phases should be even smaller. As has already been pointed 
out, this is in a sharp contrast to the conclusion of the elastic approximation which predicts 
that the phase differences are small but the phases themselves are close to 90°. Though some 
of the existing literature reach the same or similar conclusions, our s-channel picture for the 
origin of the almost real decay amplitudes is quite different from and orthogonal to them. 
Our phase-amplitude dispersion relation corroborates the smallness of the decay phases. Up 
to the short-distance QCD corrections, the magnitude of the final-state interaction phase 
is independent of isospins in the first-order approximation. This small and common final- 
state interaction phase is the reason why the color suppression appears to hold well in the 
S-decay. 

As for a phase difference between a pair of decay amplitudes, the portion of the relative 
phase that arises from the first-order correction to the random phase limit is not calculable, 
but small when the channel number n is large. The typical branching fraction for the two- 
body modes B — >• DM, D M ■ ■ ■, where M is vr, p, ■ ■ -, is somewhere between 0.1% and 1% 
while the inclusive nonleptonic branching due to b cud should be close to 70%. Therefore 
the effective number of open channels n with this set of quantum numbers is a hundred or 
more. For the decay modes such as B ^ D~p^ that have relatively large branching fractions, 
the randomness approximation gives reliable predictions while the phases of the suppressed 
decay modes such as B^ D cannot be reliably predicted in this approximation. In 
other words, the phase of the B^ ^ D amplitude can be large. The reason is that the 
two isospin amplitudes (-Dp)i/2,3/2 oi B^ ^ D p° nearly cancel each other in the B^ ^ D p^ 
combination, thus enhancing the phase fluctuation contributions. 

The relative phase which arises from the interference with the subdominant nondiffractive 
eigenchannels is negligibly small since the contribution of the nondiffractive channels is 
suppressed by /s/m^. According to our estimate in Eq. (|6TD , the phase difference of this 
origin is expected to be 

A6 ~ 1/300. (71) 

The phase difference due to the fluctuations is more important in comparison. We are unable 
to estimate the relative importance between the fluctuations and the short- distance effects 
to the phase differences. Notable exceptions are the decay modes for which the annihilation 
or exchange process dominates. For example, the decay B^ — > K^K~ for which no spectator 
diagram can be drawn. In this case the decay proceeds only through 5 — > i? — > K^K^ in 
the picture of the dual resonance model. There are no diffractive eigenchannels in the s- 
channel for this mode. Therefore the B^ K^K~ amplitude can have a large phase unlike 
the amplitude for the spectator-dominated mode B^ K^K . The rescattering from K^K 
to the K^K~ final state can occur through a qq pair exchange, namely, the non-Pomeron 
exchange. By stretching the corresponding quark diagram out, however, we find that the this 
rescattering decay process is actually due to one of the annihilation diagrams. Unfortunately 
the branching fraction is too small for such interesting decay modes. 
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VI. CONCLUDING REMARKS 



A reliable computation of the long-distance final-state interaction phases is nearly an 
impossible task. The elastic rescattering approximation is probably not viable even for 
the D decay. It is out of question for the B decay. Nonetheless we are able to extract a 
few relevant pieces of information with the phase-magnitude dispersion relation and with 
the eigenphase analysis. The conclusion from the dispersion relation is rigorous while the 
argument based on the eigenphase shifts resorts to the random channel-mixing postulate 
and the dual resonance model. The complexity of long-distance strong interactions appears 
formidable if we try to go any step further along this line. It may be possible to reach 
essentially the same qualitative conclusions with more an intuitive qualitative reasoning 
P7| . We have tried in a way different from anybody else, maybe in a hard way. We hope 
that our approach sheds a light on some aspects of the problem which have so far not been 
appreciated. 
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FIGURES 




(a) (b) 



FIG. 1. Duality in the quark diagrams for the elastic boson-boson scattering, (a) The resonant 
channels dual to the non-Pomeron Regge trajectories. Being made of a qq pair, the intermediate 
resonant states cannot have the exotic (non-gg) quantum numbers, (b) The nonresonant channels 
dual to the Pomeron consist of qqqq. In the Pomeron exchange process, one boson (gg) exchanges 
gluons with the other boson (gg). Since gluons carry no flavors, the Pomeron is necessarily a flavor 
singlet. The s-channel intermediate states consist of a pair of hadrons or hadron resonances which 
does not resonante. 
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